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ABSTRACT 

An observational model is presented for the spectra of young rotation-powered pulsars and their 
nebulae based on a study of nine bright Crab-like pulsar systems observed with the Chandra X-ray 
observatory. A significant correlation is discovered between the X-ray spectra of these pulsars (PSRs) 
and that of their associated pulsar wind nebulae (P WNe) , both of which are observed to be a function of 
the spin down energy, E. The 2 — 10 keV spectra of these objects are well characterized by an absorbed 
power-law model with photon indices, T, in the range of 0.6 < TpsR < 2.1 and 1.3 < Tpwn < 2.3, for the 
pulsars and their nebulae, respectively. A linear regression fit relating these two sets of indexes yields, 
T PWN = 0.86 ± 0.20 + (0.72 ± 0.13) x r PSR , with a correlation coefficient of r = 0.96. The spectra of 
these pulsars are found to steepen as T = T max + aE^ 1 / 2 , with T max providing an observational limit 
on the spectral slopes of young rotation-powered pulsars. These results reveal basic properties of young 
pulsar systems, allow new observational constraints on models of pulsar wind emission, and provide a 
means of predicting the energetics of such systems when lacking detected pulsations. 

Subject headings: pulsars: general — supernova remnants — stars: neutron — X-rays: general — 
radiation mechanisms: general — acceleration of particles 



1. INTRODUCTION 

The classical picture of a pulsar wind nebula (PWN) is 
that of a bright, centrally condensed, diffuse nebula whose 
broad-band non-thermal continuum spectrum is character- 
ized by a power-law with one or more breaks (see Chevalier 
1998 and Arons 2002 for reviews of PWN models). Recent 
observations of these objects with the Chandra X-ray ob- 
servatory has resolved out complex structure involving co- 
aligned toroidal arcs, axial jets, and wisps on arc-second 
scales (e.g., Crab Nebula, Weisskopf et al. 2000; Vela Neb- 
ula, Helfand et al. 2000; 2001). This basic morphology 
appears to be common among young, energetic pulsars 
associated with supernova remnants (Gotthclf 2001). Fur- 
thermore, Chandra imaging-spectroscopy in the 2 — 10 keV 
band reveals a distinct spectral signature between these 
structures and that of the central pulsar, blended in ear- 
lier studies (Gotthelf & Olbert 2002). It is now clear that 
current models for the classic PWN are inadequate in de- 
scribing their observed structure or spectrum. 

This Letter presents a spectral analysis of nine pulsars 
observed by Chandra which show bright, resolved PWNe. 
The high resolution Chandra data allows for the first time 
a systematic X-ray spectral study of pulsars and their neb- 
ulae independently. For the 2 — 10 keV energy band, the 
spectral slopes of these pulsars and their associated neb- 
ulae are found to be tightly correlated over the full range 
of measured values implying a fundamental spectral re- 
lationship not easily explained in the framework of cur- 
rent models of pulsar emission. The spectral slopes are 
also shown to follow an inverse square-root relation with 
respect to the spin-down energy - the spectral slope of 
younger, more energetic pulsars tend to be steeper, sug- 
gesting that the physics of a pulsar's particle acceleration 
is governed primarily by its spin-down energy. 

Collectively, these results provide new observational in- 



sight into the emission mechanisms of young rotation- 
powered neutron stars, characterize a pulsar's spin-down 
evolution, and strongly constrains pulsar emission models. 

2. NEW CHANDRA OBSERVATIONS 

The Chandra X-ray Observatory (Weisskopf, O'Dell, 
& van Speybroeck 1996) has targeted most pulsars with 
known X-ray bright wind nebulae. Table 1 presents the 
sample used in this study. Each pulsar was imaged by 
the Advanced CCD Imaging Spectrometer (ACIS) whose 
arc-second resolution allows us to isolate the pulsar emis- 
sion from that of the nebula. ACIS is sensitive to X-rays 
in the 0.2-10 keV spectral band with an energy resolu- 
tion of AE/E ~ 0.1 at 1 keV. All data were corrected 
for CTI effects by reprocessing the Level 1 data and ap- 
plying the standard Level 2 filtering criteria following the 
Townslcy method (Townsley et al. 2000). Time inter- 
vals of anomalous background rates associated with par- 
ticle flare events were further rejected. Data reduction 
and analysis were accomplished using the CIAO software 
package (CIAO 2.2/CALDB v.2.9). 

For each object, we extracted spectra from the pulsar, 
nebula, and their local backgrounds, when available, or 
obtained spectral parameters from the literature as noted 
in Table 1. Pulsar extraction regions are based on the im- 
age above 4.0 keV, where the count rate per pixel from the 
pulsar dominates over the surrounding diffuse nebula emis- 
sion. Large variations in nebula size and background rates 
for the sample precluded the use of a standard extraction 
aperture. Instead, the nebula region was defined by the 3er 
contour above the background level with the pulsar region 
excluded. All spectra were grouped to a minimum of 50 
counts per spectral bin and model spectra fitted using the 
XSPEC spectral fitting package (version vll.l). The ACIS 
instrument response matrix provided with the Townsley 
CTI correction software were used in these fits, along with 
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an ancillary telescope response function created accord- 
ing to standard CIAO 2.2 procedures, using the Townsley 
QEU calibration file. 

An important consideration in our analysis is spectral 
distortion due to CCD photon pile-up. Pile-up occurs 
when more than one photon is recorded as a single pho- 
ton event in a CCD read-out frame. This can severely 
skew the spectrum by registering the summed energy of 
two photons as one; it also causes the loss of soft photons 
when the recorded energy of two or more photons exceed 
the valid range of the instrument. For our measured pul- 
sars, the expected magnitude of pile-up spans the range 
of 3 — 8%. Pile-up can result in the spectral flattening of 
the pulsar spectra relative to that of the nebulae and is 
accounted for using the pile-up model available in XSPEC 
on all pulsar spectral fits. 

The extracted spectra were fitted with an absorbed 
power-law model whose flux is given by F(E) oc 
e -cr(E)Nn £j-r ^ wcrc p j s t ne photon index, Ah the inter- 
stellar hydrogen column density to the source, and <r(E) 
the Wisconsin interstellar absorption cross-section (Morri- 
son & McCammon 1983). In order to avoid potential soft 
thermal emission contamination, the spectral fits for both 
the pulsars and nebulae were restricted to energies above 
2 keV, where the fits are largely insensitive to Ah- In these 
fits, the Ah was held fixed to the value determined from 
fits to the nebula spectra over the full ACIS band. The 
spectral fits to the pulsar emission included the addition 
of a CCD pile-up model. 

Table 1 lists the resultant spectral slopes for each ob- 
ject, including the spectra slope of the pulsed emission only, 
taken from published values, when available, or derived by 
phase-resolved spectroscopy using ASCA or RXTE data. 
A comparison between pulsed and unpulsed emission mea- 
surements is a good test for systematic spectral distortion, 
as the former measurements, based on phase-resolved spec- 
troscopy using non-CCD detectors, are immune to pile-up. 
Although the spectra of the pulsed and total emission need 
not be the same (e.g., Crab; see Pravdo et al. 1997), it is 
reassuring that most agree to within measurement errors 
(see Figure 1, Gotthelf 2003). Interesting new exceptions 
are the N157B and 3C58 pulsars which predict a phase 
dependence in their spectral slopes. 

3. RESULTS 

A comparison of the spectral indices and three interest- 
ing parameters among the objects shown in Table 1 reveals 
some remarkable trends. A plot of the index for each pul- 
sar versus that of its structured nebula is shown in Figure 
1. A linear regression fit to this data, taking into account 
the uncertainties in both coordinates, yields, 

r PWN = 0.86(0.20) + 0.72(0.13) x r PSR (1) 

where the one-sigma errors are given in parentheses. This 
relationship is highly significant, with a linear correlation 
coefficient of r = 0.96, and covers the range of known pul- 
sar spectral power-law indexes in the 2—10 keV X-ray 
band. This result is rather surprising, since no such cor- 
relation was previously known or predicted. 

Another unexpected result is found by comparing the 
rank-ordered pulsar indices with the spin-down energies 
of the associated object (see Table 1). Allowing for the 



uncertainties in Tpsr, an unambiguous trend is evident, 
with a clear spectral steepening with spin-down energy 
E. Since the spin-down energy is expected to play an im- 
portant role in the evolution of the pulsar and its nebula, 
various functional forms involving the spin-down energy 
were examined in order to model the spectral slope. The 
best fit is obtained with an inverse square root model, 
T = T max + aE^ 1 / 2 , with the following parameters (see 
Figure 2): 

r PSR = 2.08(0.07) - 0.029(0.003) E m 1/2 (2) 

where E40 is the spin-down energy in units of 10 40 erg s _1 
and the one-sigma errors are given in parentheses. The fit 
is rather poor with \ 2 = 17 for 7 DoF, but nearly all the 
excess contribution to x 2 is from the Kes 75 data point, 
whose E has been previously noted to be at odds with the 
properties of the SNR (Helfand, Collins, Gotthelf 2002). 
Without this data point, the fit is excellent (x 2 — 5 for 6 
DoF), however the model parameters are not significantly 
changed and we quote the fit using all nine data points. 

Given Equation 1, we expect r P wN to follow a similar 
functional form. This is found to be the case but with a 
lesser significance of r = 0.90 (x 2 = 54 for 7 DoF), again 
with the Kes 75 point providing the greatest contribution 
to x 2 - We thus use Equations 1 and 2 to derive the comple- 
mentary equation for the spectral slopes of the structured 
nebula: 

T PWN = 2.36(0.33) - 0.021(0.005) £ 40 1/2 (3) 

The above two equations, along with the updated Se- 
ward & Wang (1988) L x /E relationship (Possenti et al. 
2002), relate the timing properties of a rotation-powered 
pulsar to its X-ray spectra in the 2 — 10 keV energy band. 
These equations predict a maximum value for the spectral 
slopes, T max , for pulsars obeying these equation. Further- 
more, in lieu of any high-energy spectral cut-off, energetic 
considerations restrict the minimum spectral slopes, and 
therefore E, to < r P sR < 2.2 and < T P wn 

< 2.7. Ad- 
ditional emission-model dependent restrictions (e.g., puta- 
tive particle distribution spectral slope; see §4) may fur- 
ther constrain the lower limit on T. 

Figure 3 shows the permitted range of r P sR in the P — P 
diagram for a young rotation-powered pulsar. Radio pul- 
sars which lie in this region, but not already in our sample, 
likely contain undetected nebulae whose spectral slopes 
can be predicted. This population may also offers a mean 
to determine the r P sR cut-off for Crab- like pulsars. The 
former is a test of the hypothesis that observable PWN 
are restricted within the P — P plane and the latter will 
provide an important constraint on pulsar emission mod- 
els. 

A trend is also noticed for the spectral hardness in the 
2 — 10 keV band which decreases with characteristic age 
r (see Table 1), but unlike for the spin-down energy, the 
correspondence is not consistently monotonic. Consider- 
ing the ordering by spin-down energy, large deviations by 
an individual pulsar from the trend may be a good indi- 
cator that t is in fact a poor measure of that pulsar's true 
age. And finally, the magnetic field is also expected to 
play an important role in the evolution of a PWN, partic- 
ular in determining the synchrotron lifetime of accelerated 
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particles (see §4). However, no simple trend is found for 
the derived magnetic field with respect to the measured 
spectral slopes or to the other calculated parameters. 

4. DISCUSSION 

The broadband spectra of the Crab pulsar and neb- 
ula have been known for some time (e.g. see Fig. 4-2 
of Manchester & Taylor 1977) and agrees well with the 
spectral relationship presented herein. But until now, no 
systematic correlation linking a pulsar's radiant spectrum 
with that of its nebula or to its spin-down energy was ei- 
ther noted or predicted. A key question is where does the 
pulsar emission occur and how is it related to the observed 
structured nebula; this is consider below in the context of 
current theory. 

Both competing models for generating high-energy emis- 
sion from a pulsar, the polar-cap (e.g., Harding et al. 1978) 
and outer-gap (e.g., Cheng, Ho, & Ruderman 1986) mod- 
els, provide a natural particle accelerator, radiating syn- 
chrotron emission close to the star (within the light cylin- 
der) in a strong magnetic field. The standard theory of 
synchrotron acceleration of particles with large pitch an- 
gles ($ pitch » 7 is the Lorentz factor) predicts a 
radiation spectrum of l v oc jy _ (p -1 )/ 2 for a particle distri- 
bution N(-y) = N -f~ p , radiating away their energy on the 
order of seconds. Most models estimate and/or simulation 
p-2 (e.g., Crusiu-Watzcl, Kunzl, & Lesch 2001). Ulti- 
mately, in a pure synchrotron model, energy constraints 
require a particle index p (= 2T — 1) greater than 2, forc- 
ing a lower bound on T > 2.0, above the synchrotron cool- 
ing frequency, or T > 1.5, below, generally in the X-ray 
regime. This is problematic for these models given the 
range of Tpsr presented in Table 1. The addition of a 
curvature radiation component and the location of the cy- 
clotron turnover may allow for spectral slope variations in 
the 2-10 keV band (e.g., Rudak & Dyks 1999). But how 
this might be connected to a pulsar's spin-down energy is 
unclear. 

In contrast to the highly localized accelerator in the 
above pulsar emission models, the standard Kennel & 
Coroniti (1984a, b) model for the Crab nebula is based on 
a symmetric, isotropic wind. In this model, a young pul- 
sar loses its rotational energy predominantly in the form 
of a highly relativistic (7 ~ 10 6 ), isotropic particle wind. 
The freely expanding wind is initially invisible as it travels 
though the surrounding self-evacuated region, but eventu- 
ally encounters the ambient medium (r s ^0.1 pc) where 
it is reverse-shocked, resulting in the redistribution of par- 
ticle energies with that of the magnetic field. The visible 
nebula is manifest as the shocked particle emitting syn- 
chrotron radiation. The above symmetric model is, how- 
ever, notably inadequate in explaining the observed tori 
and jet-like structure. 

In the context of this work, these models do not pro- 
vide a physical mechanism to link the pulsar and nebula 
spectrum. Our results suggest that the structured nebula 
is either powered directly by the pulsar or both are re- 
sponding to a common mediating process. Although the 
particles associated with the pulsed emission are available 
for the shocked wind, they are not expected to retain their 
original distribution and track the pulsar's spectrum. In 
fact, the Kennel & Coroniti (1984a, b) model for the Crab 



is inconsistent with a particle wind originating from the 
pulsar based on energetic grounds. A wind with a time- 
average injection rate of N± ~ 3 x 10 40 for pairs with 
7 ~ 10 6 would exceed the spin-down power of the Crab 
by several orders of magnitude (Gallant et al. 2001). Fur- 
thermore, in the 2 — 10 keV band, Fermi acceleration of 
highly relativistic particles yield injected spectrum with 
fixed p ~ 2.2 - 2.3 (Achterbers et al. 2001). Both the 
range of observed PWN slopes and their correlation with 
the PSR emission strongly constrain all models involving 
shock acceleration as the origin of the particle distribution. 
Investigations of collisionless models may be fruitful (e.g. 
Blanford 2003). 

The apparent steepening of the spectral slope with spin- 
down energy suggests that the spectral evolution is best 
parameterized by E instead of the characteristic age r. Of 
all the derived quantities, E is the least model dependent, 
likely to be an accurate measure of the current energy loss 
rate, r is evidentially an uncertain estimate of a pulsar's 
true age and it is less well correlated with spectral slope. 
Young pulsars are known to undergo episodes of glitches, 
and both internal and external torques can adjust their pe- 
riod and period derivatives. This allows for the possibility 
of distinct episodes of increased efficiency in the E energy 
loss during a pulsars' early evolution. This is might be 
consistent with the relatively large characteristic age for 
the bright, diffuse nebulae of N157b and 3C58, perhaps 
currently in transition. 

The PWN of some pulsars show a gradual spectral steep- 
ening away from the center of the nebulae, most notably 
for the Crab and 3C58 nebula. This is generally attributed 
to synchrotron losses from a centrally injected power-law 
distribution particle wind. If the synchrotron loss time is 
less than the pulsar's age, a steepening of Ar = 0.5 in the 
spectrum is expected across the nebula. Observationally, 
Equation 1 rules out a simple synchrotron cooling mecha- 
nism (i.e., Ar = r PWN - r PS R = -0.028 r PWN + 0.86). 
However, because of the relatively large uncertainties on 
the spectral indices, a Ar = 0.5 solution is not formally 
excluded, as the fit for the latter model is only marginally 
worse then for the former (\ 2 = 8 vs. 5, for 7 DoF, re- 
spectively). Deeper Chandra observations, with reduced 
measurement errors, are needed to distinguish fits. 

The well defined spectral relationship for our pulsar 
sample suggests that the observational properties of Vela 
as compared to the Crab are simply the consequence of 
an older, less energetic system. It is also apparent that 
the manifestation of a shell remnant associated with a 
PSR/PWN system, the main observational difference be- 
tween Vela and the Crab, must be to some degree au- 
tonomous of either pulsar or SNR evolution. This raises 
questions for the importance of a pressure confined wind, 
inherent in PWN evolution models, to account for the ob- 
served morphology of the structured nebula. 

Considering that Chandra is the first telescope to spa- 
tially resolve neutron stars from their structured nebulae, 
previous spectral observations of young pulsars are likely 
to be contaminated by nebula emission. For example, we 
see no clear relationship between spectral indices of pul- 
sars obtained by ASCA with those measured by Chandra. 
Fits to these curved composite spectra using single ab- 
sorbed power-law models, typically assumes, can result in 
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exaggerated absorption column measurements. 

For pulsars lacking measured pulsations, the relation- 
ships presented herein provide a powerful tool for predict- 
ing the spin-down energy of pulsars directly from their 
spectra. For an unresolved pulsar (and nebula) the pres- 
ence of a nebula can also be predicted, based on the mea- 
sured spectrum. As more young rotation-powered pulsars 
are observed with Chandra, a more complete picture will 
emerge of the morphological and spectral characteristics of 



young pulsar evolution, paving the way for future, perhaps 
unified, models of pulsars and their structured nebulae. 
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X-ray Spectra of Young Pulsars with Bright Central Nebulae Observed with Chandra a 



Remnant 


Pulsar 


1 PWN 


pc 

1 PSR 


1 Pulsed 


(kyr) 


E e 

(xlO 35 ergs/s) 


Bp/^QED 


Gil. 2-0.3 


PSR J1811-1926 


1.28±0.15 


0.63±0.12 


0.60±0.60 


23.0 


53 


0.04 


Vela XYZ 


PSR J0835-4510 


1.50±0.04 


0.95±0.24 


0.93±0.26 


12.0 


67 


0.08 


G54.1+0.3 


PSR J1930+1852 


1.64±0.18 


1.09±0.09 


1.06±0.86 


2.9 


118 


0.23 


Kes 75 


PSR J1846-0258 


1.92±0.04 


1.39±0.11 


1.10±0.30 


0.7 


82 


1.10 


MSH 15-52 


PSR J1513-5908 


1.93±0.03 


1.40±0.50 


1.26±0.08 


1.6 


140 


0.34 


3C 58 


PSR J0205+6449 


1.92±0.11 


1.73±0.15 


1.11±0.34 


5.0 


263 


0.08 


SNR 0540-69 


PSR J0540-6919 


2.09±0.11 


1.88±0.11 


1.83±0.13 


1.7 


1481 


0.11 


Crab Nebula 


PSR J0534+2200 


2.14±0.01 


1.85±0.09 


1.87±0.05 


1.3 


4394 


0.08 


N157B Neb. 


PSR J0537-6910 


2.28±0.12 


2.07±0.21 


1.60±0.35 


5.0 


4916 


0.02 



a Ranked by increasing pulsar photon index, TpwN (then Tpgj^, were ambiguous). Data for the Crab pulsar and 
nebula were obtained from Pravdo, Angelini & Harding (1997). 

b Values for the following objects were taken from the literature: G54.1+0.3: Lu et al. 2002; SNR 0540-69: Kaaret 
ct al. (2001). 

c Includes both pulsed and un-pulsed emission, corrected for pile-up. 

d Pulsed spectrum references: Vela: Strickman, Harding & Jager (1999); Gil. 2-0. 3: Torii et al. (1997); Kes 75: 
Gotthelf et al. (2000); N157B: Marshall et al. (1998); SNR 0540-69: Kaaret et al. (2001). 

e The characteristic pulsar spin-down age is defined as r = P/2P, the spin-down energy as, E = 4,n 2 IP/P 3 , were 
/ = 10 45 gm cm~ 2 , and the inferred magnetic field B p = — 6c 3 _B/_R 6 f! 4 G is normalized to the quantum critical field, 
B QED = mlc z /eh = 4.4 x 10 13 G. 
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Fig. 1. — Relationship between a pulsar's spectral slope (Tps^) and that of its structured nebulae (rpwr-j) m the 2 — 10 keV energy range, 
assuming a simple power-law model for the objects presented in Table 1. The dashed-line indicates the best-fit linear model. The physical 
origin of this relationship has yet to be determined 
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(E/10 4U erg/s) 

Fig. 2. — A comparison between the spectral slope of the 2 — 10 keV pulsar emission (Fpsr) and the square root of the spin-down energy, 



E] 



-1/2 



in units of 10 40 erg s 1 , for the each object presented in Table 1. The dashed-line indicates the best-fit model. 
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PULSAR PERIOD (SEC) 

Fig. 3. — The period-period derivative diagram showing the sample of X-ray objects presented in Table 1 (open circles), a complete sample 
of known radio pulsars from Manchester (2003) (filled circles), and the known AXPs/SGRs (stars). Lines of constant spectral slope for a 
young, rotation-powered pulsar (in the 2 — 10 kcV energy band; see text) is given by the diagonal lines. The solid line corresponds to Tpsr = 0. 



